Fidelity Decay as an Efficient Indicator of Quantum Chaos 
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Recent work has connected the type of fidehty decay in perturbed quantum models to the presence 
of chaos in the associated classical models. We demonstrate that a system's rate of fidelity decay 
under repeated perturbations may be measured efficiently on a quantum information processor, and 
analyze the conditions under which this indicator is a reliable probe of quantum chaos and related 
statistical properties of the unperturbed system. The type and rate of the decay are not dependent 
on the eigenvalue statistics of the unperturbed system, but depend on the system's eigenvector 
statistics in the eigenbasis of the perturbation operator. For random eigenvector statistics the decay 
is exponential with a rate fixed precisely by the variance of the perturbation's energy spectrum. 
Hence, even classically regular models can exhibit an exponential fidelity decay under generic quan- 
tum perturbations. These results clarify which perturbations can distinguish classically regular and 
chaotic quantum systems. 

PACS numbers: 05.45.Mt, 03.67.Lx 



Over the last two decades a great deal of insight has 
been achieved regarding the manifestations of chaos and 
complexity in quantum systems. We are interested in the 
problem of identifying such signatures in the context of 
quantum simulation on a quantum information processor 
(QIP). It is known that a QIP enables efficient simulation 
of the dynamics of a wide class of quantum systems 
In the case of quantum chaos models, quantum simula- 
tion provides an exponential speedup over direct classi- 
cal simulation 4, 5]. Recently, the quantum baker's 
map has been implemented using a nuclear magnetic res- 
onance QIP . These developments highlight the impor- 
tance of devising efficient QIP methods for the measure- 
ment of quantum chaos signatures and related properties 
of complex quantum systems. 

Perhaps the most established signature of quantum 
chaos is given by the (nearly) universal correspondence 
between the eigenvalue and eigenvector Q statistics 
of quantized classically chaotic systems and those of the 
canonical ensembles of random matrix theory (RMT). 
Unfortunately direct detection of these spectral signa- 
tures is algorithmically inefficient by any known tech- 
nique. However, following the original observation of 
Peres 0, some recent work has demonstrated that, un- 
der sufficiently strong perturbation, the fidelity exhibits 
a characteristic exponential inthe case of classical chaotic 
systems 0, [ul 111, [S] ■ Below, we demonstrate that this 
charecteristic decay may be measured by an efficient algo- 
rithim, and analyze in detail how the observation of this 
decay may be applied as an indicator of canonical RMT 
statistics (quantum chaos) in the unperturbed system. 

By first considering RMT models we show that 
Wigner-Dyson fluctuations in the system eigenvalue spec- 
trum are not necessary to produce this characteristic 
decay. More importantly, we find that the canonical 



RMT statistics for the system eigenvector components, 
expressed in the eigenbasis of the perturbation operator, 
are sufficient to produce the characteristic exponential 
decay. These observations are checked in the case of a 
dynamical model, where we demonstrate that an expo- 
nential fidelity decay can arise not only for a classically 
chaotic system, but also for classically regular system un- 
der a generic choice of perturbation. Our analysis implies 
specific restrictions that must be imposed on the choice of 
applied perturbation operator in order to extract useful 
information about the statistical signatures of quantum 
chaos in the unperturbed system. 

The aim of this Letter is to characterize certain static 
properties of a unitary map, U , by observing the rate 
of divergence between a fiducial Hilbert space vector 
evolved under this map, \ipu{n)) = C/"|'(/'(0)), and the 
same initial vector evolved under this map but sub- 
ject also to a sequence of small perturbations \ilip{n)) — 
(t/pC/)"|V'(0)), where Up = exp(— iV^(5) is some unspeci- 
fied perturbation operator and n denotes the number of 
iterations. The fidelity. 



0(n) = |(V'„(n)|Vp(n))|^ 



(1) 



provides a natural indicator of this divergence. The 
value of 0{n) may be determined by an efficient al- 
gorithm on a QIP as follows. We start by preparing 
the fiducial state |V'(0)) = C/qIO)®"", where is the 
number of qubits required to span the system's Hilbert 
space N — 2"'. As shown below, the choice of fidu- 
cial state is not critical and the computational basis 
states, which are simplest to implement, are a conve- 
nient set. After applying the sequence ([/'^ )"([/pJ7)", the 
system register contains a final state \'ip{n)). The circuit 
implementation of U requires only Poly(rig) operations 
for the simulation of a wide class of quantum systems 
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to arbitrary accuracy 0, S S H Here we observe 
that KVpHIV-^n))! = mnMiO))\ = \{^f\0)l where 
the state \tjjf) is obtained by time-reversing the initial 
state preparation \ipf) = U^\ip{n)). The magnitude of 
0(n) is then determined from samphng the population 
of the state |0). The entire algorithm therefore scales as 
Poly (rig). 

Recently, Jacquod and coworkers [llj observed that 
for I "0(0)) = I Wo), an initial eigenstate of U with eigen- 
phase (j)o, the fidelity relates to the local density of states 
(LDOS), r]{(j)o — = |(wo|Wm)P, via Fourier transform, 



0{n) = I V?7((/)o -(/)'„) exp(-i((/io-(/)J„)n) 



(2) 



In the above jw^) are eigenstates of the perturbed map 
UpUlv'^) = exp{—i(j)'^)\v'^) . These observations locate 
the decay of 0(n) in an existing theoretical framework. 
In the non-perturbative regime cr/A > 1, where cr^ = 
S'^V^^ denotes a typical off-diagonal matrix element and 
A is the average level spacing, previous work suggests 
that when the perturbed system is complex the LDOS is 
typically Lorentzian [ij, [l^ , 



0JJ^ + (r/2)^ 



(3) 



with width F = 27r cr^ /A determined by the Fermi golden 
rule (FGR) . From (|2J) and Q one expects the exponential 
decay. 



0{n) ~ exp(— Fn). 



(4) 



The onset of the exponential decay I0J has been con- 
firmed recently in a few classically chaotic models 0, 
ll^ , though the rate is not always given by the golden 
rule fl^ . However, in the case of integrable U the 
situation is less clear since under some perturbations the 
LDOS is known to take on a Lorentzian shape jlfil. 11^. 

Below, we examine which statistical properties of the 
unperturbed system lead to the FGR decay, and how this 
characteristic decay depends on the properties of the per- 
turbation operator. This approach is motivated by the 
context of QIP simulation, in which the eigenbasis of this 
perturbation may be mapped onto an arbitrary basis of 
the simulated system U • Our choice of perturbation 
eigenvalue structure is motivated from the perspective of 
quantum control studies. Specifically, we consider 



C/p = n"^iexp(-ifoi/2), 



(5) 



where az is the usual Pauli matrix and Q therefore cor- 
responds to a collective rotation of all the qubits b y an 
angle S. Eq. |3l is a model of coherent far-field errors 
and for this type of error model a better understanding 
of the fidelity decay is a subject of intrinsic interest. 



As a first test of the LDOS/FGR framework in the case 
of the qubit perturbation (0) we evaluate the fidelity de- 
cay for a map U = J7cue drawn from the circular unitary 
ensemble (CUE). These matrices form an established 
model for classically fully chaotic time-periodic systems 
(without additional symmetries) since these systems (al- 
most always) exhibit the same characteristic (Wigner- 
Dyson) spectral fluctuations and eigenvector statistics as 
the CUE [i^. The randomness of the system eigenvec- 
tors enables a system-independent estimate of the rate 
F of the FGR decay. Since the components of random 
eigenstates are distributed uniformly over the basis states 
and uncorrelated with the distribution of eigenvalues, the 
second moment of the matrix elements Vmn may be di- 
rectly evaluated. 



(6) 



assuming A = and where — X]t=i denotes 
the variance of the eigenvalues of V. As a result, the rate 
of the FGR decay is determined by the eigenvalues of the 
perturbation, 



A2, 



(7) 



where we have used A = 2Tr/N. For the qubit pertur- 
bation ((nj the variance of the eigenvalues has a simple 
form. 



A2 



N ^ 

k=0 



2k — rir; 



(8) 



where the C"' are binomial coefficients. Using our RMT 
estimate O, for Uq — 10, the rate is. 



F = 2.50 



(9) 



While a CUE map may be generated on a QIP using 
the gate decomposition devised in Ref. [U, for our nu- 
merical study we construct U = C/cue directly from the 
eigenvectors of a random Hermitian matrix. Since com- 
putational basis states are easiest to prepare in the QIP 
setting, we consider the fidelity decay for both single com- 
putational basis states and averages over 50 such states. 
The behaviour of the fidelity decay for a matrix typical 
of CUE is displayed in Fig. ^ The three perturbation 
values displayed in the figure are chosen near the onset 
of the non-perturbative regime (8 > 0.1) and it is evident 
that the fidelity decay even for individual computational 
basis states exhibits FGR decay (0J) at the expected rate. 
The FGR decay persists for a time-scale F~^ log(iV) until 
saturation at a time-average that decreases as 1/A^. 

We next consider 0{n) for the Gaussian unitary ensem- 
ble (CUE) in order to clarify the relationship between 
the FGR decay and the distinct statistical features of 
RMT that represent signatures of quantum chaos. The 
CUE consists of Hermitian matrices with independent 
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FIG. 1: Fidelity decay for t/cuE averaged over 50 compu- 
tational basis states (dash lines) is in excellent agreement 
with the golden rule decay (|1J and the RMT rate @ for 
S — (0.1,0.2,0.4) (solid lines top to bottom). Chain lines 
show the fidelity decay for two typical computational basis 
states in the case S — 0.1. 

elements drawn randomly with respect to the unique 
unitarily-invariant measure [2i| . GUE forms the relevant 
RMT model for the important class of chaotic or complex 
autonomous Hamiltonian systems that are unrestricted 
by any additional symmetries. We may examine the sen- 
sitivity to perturbations for the GUE by constructing 
the unitary operator C/gue — exp(— ii/GUET"), where r 
is a time-delay between perturbations. We consider the 
same perturbation as for the CUE case. For sufficiently 
small T the propagator approaches identity and the over- 
lap decay is dominated by the perturbation operator, 
0{n) = |(exp(-i(5y))|2 + ©((JV^n^). This behaviour is 
demonstrated in Fig. |21for r — 0.001 and r = 0.01 and 
with S = 0.3. For larger values of r the fidelity decay 
under Ugve obeys the FGR with the RMT rate ©. 

The important point is that for sufficiently large t the 
eigenphases of J7gue become spread pseudo-randomly in 
the interval [0,27r). Under these conditions the eigen- 
phases of the map C/gue exhibit the Poissonian spec- 
tral fluctuations that are characteristic of classically inte- 
grable (time-periodic) systems. We checked the nearest- 
neighbor spacing distribution of C/gue and found that 
for T = 100 the statistics are in excellent agreement with 
the Poissonian distribution P{s) oc exp(— s) (see inset 
to Fig. 12)). However, the eigenvectors of C/gue are ran- 
dom (by construction) and independent of t (for finite r) . 
From these observations it is clear that the presence of 
Wigner-Dyson spectral fluctuations in the implemented 
U , which comprises the only basis- independent criterion 
of quantum chaos, is not actually necessary for the onset 
of exponential (FGR) decay at the rate Q. This sug- 
gests that it is the RMT statistics of the eigenvectors of 
U that lead to the FGR decay with the RMT rate. 

We next consider the fidelity decay for the quantum 
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FIG. 2: Fidelity decay for C/gue averaged over 50 computa- 
tional basis states with S = 0.3 for r — 0.001 (dashed line), 
r — 0.01 (dotted line) and r = 0.1 and 100 (chain lines), com- 
pared to the FGR/RMT prediction (solid line). Inset: Spac- 
ing distribution for C/gue for r = 100 compared to Poissonian 
(solid line) and Wigner-Dyson (dashed line) distributions. 

kicked top, which is an exemplary dynamical model of 
quantum chaos [20I E^ . The kicked top is a unitary 
map C/qkt = exp(— z7rJy/2) exp(— zfcJ|/j) acting on the 
Hilbert space of dimension N = 2j + 1 associated with 
an irreducible representation of the angular momentum 
operator J. In previous fidelity decay and LDOS stud- 
ies the choice of perturbation has usually been tied to a 
physical coordinate of the system U. We first follow this 
convention and identify the eigenbasis of the perturbation 
(O with the eigenbasis \mj) of the system coordinate Jz 
(where rrij = {j, . . . , —j}). In Fig. Owe compare the fi- 
delity decay for the chaotic and regular regimes of the 
kicked top for averages over 50 initial computational ba- 
sis states. The fidelity decay for the chaotic top {k = 12) 
is well described by the FGR prediction (gj and the RMT 
rate whereas the regular top {k — 1) shows a slower 
non- exponential decay rate. Similarly, if we associate the 
perturbation eigenbasis with the basis of the Jy coordi- 
nate, the fidelity decay for the chaotic top remains in 
agreement with the RMT rate and the regular top again 
exhibits non-exponential decay, though in this case with 
a faster decay than the RMT rate ©. However, we now 
demonstrate that an exponential decay at the RMT rate 
arises even for the regular kicked top when the qubit per- 
turbation jSJ is diagonal in a generic basis relative to the 
eigenbasis of C/qkt- Specifically, we leave the perturba- 
tion eigenvalue spectrum unchanged but set 

Up^T [Hj exp{-i5al/2)] T-\ (10) 

where T is drawn from CUE. As demonstrated in the 
inset to Fig.|3 under this type of perturbation the fidelity 
decay for the regular top is indistinguishable from that 
of the chaotic top and is very accurately described by the 
FGR at the RMT rate ®. 
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FIG. 3: Decay of 0{n) for the kicked top in chaotic regime 
{k = 12) averaged over 50 computational basis states, with 
the perturbation eigenbasis mapped to the eigenbases of 
(dash lines) and Jy (chain lines) compared to the FGR decay 
(solid lines) at the RMT rate @ for 5 = (0.1,0.3) (top to 
bottom). Lines with circles and squares are for the regular 
kicked top (fc = 1) with the perturbation eigenbasis tied to the 
Jz and Jy coordinate bases respectively (for 5 = 0.1). Inset: 
Average fidelity decay for regular kicked top (fe = 1), when 
the qubit perturbation is in a random eigenbasis 1101 and with 
S = (0.1, 0.3) (dashed fines), compared to the FGR/RMT rate 
(solid lines). 

The sensitive dependence of the type of fidelity decay 
on the eigenbasis of the applied perturbation suggests a 
close connection with the basis-dependence of the eigen- 
vector statistics of classically regular quantum models. 
Expressed in a generic quantum basis, the eigenvectors 
of any quantized classical system U will have randomly 
(Gaussian) distributed components. In contrast, in the 
eigenbases of the system coordinates the components of 
classically chaotic and integrable systems are known to 
be different with the former Gaussian distributed 
and the latter exhibiting substantial deviation from the 
canonical Gaussian distribution. In light of this connec- 
tion, in the case of quantized classical models we infer 
that exponential (non-exponential) fidelity decay can be 
correlated with the presence (absence) of characteristic 
RMT spectral fluctuations in the unperturbed system 
provided that the applied perturbation commutes with 
a system coordinate. 

In summary, we have shown that the fidelity decay 
may be measured efficiently on a QIP. We then exam- 
ined which statistical properties of the unperturbed sys- 
tem determine the type and rate of the decay. In the case 
of random unitary and Hermitian matrices, as well as a 
classically chaotic dynamical model, we have shown that 
the fidelity decays exponentially with a characteristic 
rate given precisely by the variance of the perturbation's 
eigenspectrum. The occurrence of the exponential decay 
is not directly dependent on the Wigner-Dyson fluctua- 
tions of the unperturbed spectrum, but does depend sen- 



sitively on the RMT statistics of the system eigenvectors 
in the eigenbasis of the applied perturbation. Hence, the 
fidelity decay for both classically regular and chaotic dy- 
namical systems is given by the FGR under all but a small 
subset of unitary perturbation operators. In the case 
of classical models, we conclude that the fidelity decay 
provides a reliable indicator of RMT statistics (quantum 
chaos) in the unperturbed system only when the applied 
perturbation is restricted to the subset of perturbations 
that commute with a classical coordinate. 
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cial support. 
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